Dynamic entanglement in oscillating molecules and potential biological implications 
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We demonstrate that entanglement can persistently recur in an oscillating two-spin molecule 
that is coupled to a hot and noisy environment, in which no static entanglement can survive. 
The system represents a non-equilibrium quantum system which, driven through the oscillatory 
motion, is prevented from reaching its (separable) thermal equilibrium state. Environmental noise, 
together with the driven motion, plays a constructive role by periodically resetting the system, even 
though it will destroy entanglement as usual. As a building block, the present simple mechanism 
supports the perspective that entanglement can exist also in systems which are exposed to a hot 
environment and to high levels of de-coherence, which we expect e.g. for biological systems. Our 
results furthermore suggest that entanglement plays a role in the heat exchange between molecular 
machines and environment. Experimental simulation of our model with trapped ions is within reach 
of the current state-of-the-art quantum technologies. 



The question, to what extent quantum mechanics plays 
a role in biology, is still far from being well-understood 
pi, ^ . It seems that classical concepts alone are in- 
sufficient for a proper understanding of certain biolog- 
ical processes, and that coherent quantum effects need 
to be taken into account. It has e.g. long been known 
that quantum tunneling plays an important role in enzy- 
matic reactions [3l [4j . Experimental evidence for quan- 
tum coherence in the photosynthetic system has recently 
been reported in [5]. The interplay between the coher- 
ent free Hamiltonian and the environment is believed to 
significantly enhance quantum transport in the Fenna- 
Matthews-Olson (FMO) protein complex QEIE]. 

Apart from these isolated instances where quantum co- 
herence seems to help, it is however not clear to what ex- 
tent biological systems exploit quantum mechanics e.g. 
to optimize their functionality (beyond the trivial fact 
that the latter determines, of course, the structure of 
bio- molecules) . Most physicists and biologists are gener- 
ally skeptical about the question whether genuine quan- 
tum features such as entanglement play a broader role in 
biology. The obvious reason for that viewpoint is that 
entanglement is very sensitive to noise and requires spe- 
cial conditions to be maintained, in particular very good 
insulation. Biological systems are anything but - they 
are wet and hot, and with extremely high levels of noise. 

An often ignored fact is, however, that biological sys- 
tems are also open driven quantum systems, operating far 
away from thermal equilibrium [2 . This opens many new 
possibilities which have not yet been carefully considered. 
Different from e.g. solid state physics, things in biology 
move. Protein functions, for example, require conforma- 
tional motion [lOl[TT]. During such motion, e.g. in the 
context of protein folding or in isomeric transitions, we 
have to consider time-dependent quantum interactions 
(capable of forming e.g. hydrogen or ionic bonds) which 
are effectively switched on and off while the molecule 



changes its shape. Since these interactions are accompa- 
nied by a substantial amount of noise (e.g. from fluctu- 
ating dipole fields from the hydration shell and the bulk 
solvent), they are usually treated classically. A proper 
understanding of protein dynamics may however require 
one to explore the capability of these time dependent in- 
teractions and whether they need to be treated quantum 
mechanically [12]. It is, for example, not clear whether 
or not entanglement is generated during these motional 
processes. A positive answer to this question might re- 
veal novel and subtle aspect of protein dynamics and bio- 
motoric processes. It would also provide a new twist to 
the study of nano-bio interfaces, which are just in their 
infancy. 

While it seems reasonable to describe the conforma- 
tional motion of a bio-molecule classically, it also car- 
ries quantum degrees of freedom, such as nuclear spins 
and electronic states. These give rise to effectively time- 
dependent quantum interactions, with their strengths de- 
pending on the molecular shape, see Fig. [l] We analyze 
the role of environmental noise and decoherence in such 
interactions. We find that entanglement can be generated 
even at room temperature and despite the presence of 
decoherence, suggesting that the underlying interactions 
should indeed be described fully quantum mechanically 
to account for more subtle processes. 




Figure 1: Conformational changes of a bio-molecule TP, in- 
duced e.g. by the interaction with some other chemical, can 
lead to a time-dependent interaction between different sites 
(blue) of the molecule. See also 
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As a paradigmatic example, we study the time- 
dependent interaction of two spins in a thermal and de- 
coherent environment. We may imagine that the spins 
are attached to some classical backbone structure whose 
shape changes in time, as drawn schematically in Fig. |2j 
For simplicity we call such an arrangement a two-spin 
molecule. We demonstrate that, if the distance between 
the spins is oscillating, cyclic generation of fresh entan- 
glement can persist, even if no static entanglement can 
survive. Environmental noise plays thereby both a de- 
structive and constructive role by effectively resetting the 
system [13]. The oscillating molecule may be viewed as 
a molecular machine that exchanges heat with the reser- 
voir. Our results then suggest that entanglement is rel- 
evant to the absorption of heat from the environment, 
which might possibly affect certain biological processes. 




Figure 2: Model of a two- spin molecule which undergoes con- 
formational changes as a function of time. Both the spin-spin 
interaction strength J and local fields B are position depen- 
dent. 

In the semi-quantal picture that we have introduced 
above, the conformational changes lead to classical mo- 
tion of quantum degrees of freedom as illustrated in Fig.|2j 
We assume that the two spins are coupled with Ising in- 
teraction and that there also exist local electric and/or 
magnetic fields, both of which are usually position de- 
pendent. Thus, the classical molecular motion induces 
an effective time-dependent Hamiltonian of the form 

HM{t) = J(t)a(i)a(2) + B{t){ai'^ + af^) (1) 

where cr^^ and ai"^ are Pauli operators of the ath spin, 
J it) is the interaction strength, and uoo{t) = 2B{t) the 
local level splitting. We emphasize that the subsequent 
results also hold for more general Hamiltonians, but for 
simplicity we concentrate here on the Ising interaction. 
The coupling of the spins to the environment will be de- 
scribed by a master equation of the form 

^^pit) = -i[HM{t),p] + Vp{i) ^ C{t)p{t) (2) 

where Vp = 2L^pI/]^ — Lj^L^p — pLj^Ljj^ describes the 
effect of the molecule-environment coupling, and are 
Lindblad-type generators. 

The effect of the environment on the motion of bio- 
molecules is complex and far from being understood. 



To demonstrate the essential physics, we first consider a 
worst-case scenario, where the environment is described 
by bosonic heat-bath, with each spin being coupled to 
an independent thermal bath of harmonic oscillators 
[14|[T5]. In the static case, it is well-known that, above a 
certain temperature, no initial entanglement can persist 
in such an environment. We will however show that, even 
under such unfavorable conditions, entanglement can be 
generated if the particles start oscillating and the system 
moves out of equilibrium. 

If the molecular oscillation is not too fast, in the sense 
that the adiabatic condition for closed systems is sat- 
isfied, the effect of the environment on the oscillating 
molecule can be described by a master equation of type 
([2|, with implicitly time-dependent Lindblad generators 
(see also Supplementary Information). As far as the 
static entanglement is concerned, at every molecular con- 
figuration (with fixed J and 5), the molecule will be 
driven towards its thermal equilibrium state pth at tem- 
perature T. In the following, we adopt the concurrence 
C{p) [16] as the measure of two-qubit entanglement of 
a state p. For a separable (non-entangled) state it van- 
ishes, while for a maximal entangled state it reaches the 
value 1, i.e. < C{p) < 1. It can be shown that, if the 
temperature is above a critical value Tc, no entanglement 
can survive in any static configuration of the molecule, 
i.e. C{pth) = 0. In the following, we will only consider 
such situations where T > Tc. 

The main question we are asking is this: Can entan- 
glement possibly build-up through the classical motion of 
the molecule? The answer is affirmative and we demon- 
strate that entanglement can indeed persistently recur in 
an oscillating molecule, even if the environment is so hot 
that the static thermal state is separable for all possible 
molecular configurations, i.e. T > max{Tc}. 

Let us first present an intuitive explanation. Con- 
sider the following simple process: until time t = 0, the 
spins are kept distant (with J 0) and the molecule 
is in the thermal equilibrium state, with the fraction pg 
of the population in the ground state (|eo(0)) |ii))- 
If the local level splitting is sufficiently large such that 
^^o(0) ^ ksT^ Pg will be relatively large compared to 
the other energy levels. The thermal state will therefore 
be close to the ground state, which is, in this case, non- 
entangled. The adiabatic molecular motion from the dis- 
tant configuration to proximity will transform the eigen- 
states of Hm{0) into those of HM{t); in particular, the 
ground state |eo(0)) \^o{t)) will become entangled as 
the coupling between the spins increases. This explains, 
qualitatively, why we may expect entanglement to build 
up in one run of a conformational change, given that 
the molecular motion is slow enough to be adiabatic, but 
at the same time faster than the thermalization process: 
Driven through the classical motion, the system is so-to- 
speak "kicked out" of the (separable) thermal equilibrium 
state, as can be seen in Fig. [Sj 
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The above analysis only suggests that entanglement 
may appear in one run of a conformational change. How- 
ever it does not explain how one can expect to see en- 
tanglement on a longer time scale, when the environment 
begins to mix the internal states as the molecule contin- 
ues to oscillate. It seems that, in the long run, it may 
(and will) disappear as usual. What we are interested in, 
however, is the persistent generation of dynamic entan- 
glement, thus an extra mechanism is necessary to refresh 
the state of the molecule by resetting it back to the ini- 
tial state. It is intriguing that this role can be played 
by environmental noise together with oscillatory motion, 
both of which naturally exist in biological systems with- 
out further need for control. 

Let us now present the numerical results which we ob- 
tained by numerically integrating equation [2j We con- 
sider the situation where the spin positions are 

27rf 

a^aW =x«(0) + (-l)"a(cos 1) (3) 

r 

where Xa{0) are the initial positions, a is the amplitude 
of oscillation, and r is the oscillation period. For the 
local fields we assume Gaussian functions of the spin 
position as B{t) = Bq — Bie~^ For the interac- 

tion between two spins we assume dipole-dipole coupling 
J{t) = Jo/d^{t) with d{t) = \xi{t) —X2{t)\. It can be seen 
from Fig. [3]that recurrent fresh entanglement appears on 
the asymptotic cycle. The thermalizing environment is 
here constructive by re-pumping the population into the 
ground state. 




Figure 3: (Color online) Left: Ground state populations^ and 
entanglement C vs. the molecular configuration characterized 
by spin-spin distance d for the bosonic heat bath with unit 
temperature T — 1 and the system-bath coupling strength 
/t: = 0.01. The blue curve is the instantaneous thermal equilib- 
rium state. The arrows indicate the evolution direction as the 
molecule oscillates. Right: Heat current Jh vs. entanglement 
C. The oscillation parameters are xi(0) = — X2(0) = —20, 
a = 5, T = 100, and Bq 1.3, Bi 2.4, a 120, Jo = 1 x 10^ 
(see Eq. ([3| and the text thereafter). We use dimensionless 
parameters for which the Planck and Boltzmann constants 
are set to unity, h = 1, kh = 1. The temperature T = 1 
means that we set the thermal energy khT as the unit and 
express the other values as multiples of the thermal energy. 
For example, B = 1 indicates that KB /khT = 1. 

For biological systems, T ^ 300i^, the thermal energy 



ki)T is about 0,025 eV. Energy scales in bio-molecules 
are typically of the order 0.01 — 0.5 eV (e.g. for hydrogen 
bonds or electronic excitation [l2l[T7]) and thus hx)Q{ff) 
can be several times larger than /C5T. The system-bath 
coupling strength n = 0.01 corresponds to a thermaliza- 
tion time of the order of ~ ps, which compares e.g. with 
the timescale of (fast) conformational changes within bio- 
molecules and with relaxation times in the FMO com- 
plex [18] . These numbers seem to be consistent with our 
conclusion that recurrent entanglement might indeed be 
found in bio-molecular processes at room temperature. 

The oscillating molecule exhibits rather distinct fea- 
tures of non-equilibrium thermodynamics e.g. the 
entropy does not always increase with time and reach 
its saturate value. The most intriguing feature is a con- 
nection between entanglement and the heat current J'h 
between the molecule and its environment, which is de- 
fined by the energy dissipated via the heat bath as [M] 

Mt)=Tr{HM{t)[Dp{t)]}. (4) 

It can be seen from Fig. [3] that, whenever entanglement 
appears, J'h{t) is always positive, i.e. the molecule tends 
to absorb heat. This provides some evidence that en- 
tanglement is related with the heat exchange between 
the molecule and its environment. Since the oscillating 
molecule is not in thermal equilibrium, one cannot adopt 
the standard definition of temperature. However, by us- 
ing the spectral temperature as defined in [20 , one finds 
that the molecule is effectively cooled down through the 
classical motion, even though the attached thermal bath 
is always at a fixed higher temperature. This interesting 
observation is consistent with our common intuition that 
entanglement appears as the system is cooled down. 

The results discussed so far have been obtained by 
modeling the noisy environment as an Ohmic bosonic 
heat bath. Clearly, this is a highly idealized model and 
in any real biological scenario we have far more com- 
plicated interactions, e.g. with the surrounding hydra- 
tion shell and the bulk solvent [121 [H] • We found simi- 
lar results also with other decoherence models, based on 
collision-type interactions of the environment with the 
system [19 . These can be described by Lindblad gen- 
erators Lg^^ = ^/jsa'^^ and L^J^^ = y^"y{l — s)cri"\ 
where 7 is the collision-induced effective relaxation rate 
and s is the mean excitation of a spin in thermal equilib- 
rium. In Fig. [4] we demonstrate the competition between 
the constructive and destructive effects of environmen- 
tal noise in such a model, which yields an optimal value 
for the oscillation period r to establish entanglement, see 
Fig. [4^. For short oscillation periods, efficient thermal- 
ization becomes more important: with growing rate 7, 
the increase in reset efficiency more than compensates 
the decrease in efficiency of entanglement generation, see 
Fig. ^p. In this regime, the net effect of environmental 
noise is constructive. The dependence of the entangle- 
ment on r and 7 seems reminiscent of what happens 
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in the very nice example of quantum "stochastic reso- 
nance", which has been described in a quantum-optical 
context [22l. That phenomenon is, however, fundamen- 
tally different from ours, as it involves a bath at zero 
temperature, and the noise is the main driving force. 
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Figure 4: (Color online) Left: The maximal value of entan- 
glement on the asymptotic cycle Cm vs. the oscillation pe- 
riod T (logarithmic coordinate) with 7 = 0.1. Right: Cm 
vs. 7 with T = 10. The other parameters are s — 0.2, 
xi(0) -X2(0) -25, a 10, and Bo = Bi = 1.2, a 120, 
Jo = 1.2 X 10^ 



Given the complexity of biological systems, how to 
characterize biological environment is far from been un- 
derstood. The effect we presented — i.e. the existence of 
persistent entanglement — is to a very large extent inde- 
pendent of the precise details of the classical movement, 
and the thermal bath. Of course, the detailed charac- 
teristics of the entanglement (how much entanglement, 
how does it vary with time, etc.) depend on the driving 
motion and the specific environment, but the very ex- 
istence of persistent entanglement is generic (extensible 
to different types of classical motions, spectral density 
of thermal bath and also to the non-Markovian environ- 
ment). A number of different models illustrating this 
generic property are presented in the Supplementary In- 
formation. 

One can test the feasibility of our model by simulat- 
ing an oscillating molecule in a noisy environment. Two 
internal levels of trapped ions can encode an effective 
two-level system. The system Hamiltonian in the form 
of Eq.Q is implementable via state-dependent optical 
dipole forces [2j. The classical oscillation, the essence 
of which is to introduce a time-dependent Hamiltonian, 
can be simulated by tuning the interaction strength and 
the transverse fields, which is achievable e.g. by chang- 
ing the amplitudes of laser beams [25l [26l [27] . Decoher- 
ence can be simulated by e.g. applying random pulses 
to induce different decoherence channels. One can also 
simulate the bosonic bath by engineering the coupling 
between ions and vacuum modes of the electromagnetic 
field through laser radiation [28j. Entanglement can fi- 
nally be detected by performing quantum state tomog- 
raphy as in [29]. Other implementations are conceivable 
e.g. using quantum dots mounted on the tips of oscillat- 
ing cantilevers [30j . 
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SUPPLEMENTARY INFORMATION 

This is supporting material for our paper. We present 
the essential steps in the derivation of the quantum mas- 
ter equation for the oscillating molecule in contact with 
bosonic heat baths, and the calculation of the concur- 
rence which quantifies the entanglement between the 
spins. We also describe the spin-gas model as an alterna- 
tive model for the environment and show that it gives rise 
to qualitative similar features of reset and recurrent en- 
tanglement. Finally, we analyze the competing features 
of environmental noise for the generation and degrada- 
tion of entanglement in such a model. 

Methods.— To account for the effect of the environ- 
ment on the oscillating two-spin molecule, we have stud- 
ied different models, two of which we discuss here. In 
the first model - the bosonic heat bath - each qubit is 
coupled to an independent thermal bath of harmonic os- 
cillators. This is a well-known decoherence model and 
the derivation of the master equation (L2) [11] follows 
standard techniques (see e.g. [1]), with the important 
difference that in our case the Hamiltonian HM(f) of the 
system is time dependent^ i.e. its instantaneous energy 
spectrum changes as the qubit s move. This makes the 
analysis much more complicated, but under certain con- 
ditions it can be still be described by an equation of type 
(L2), albeit with time-dependent Lindblad operators. 

In the second model - the spin-gas model [2 - each 
qubit is subject to random, collision-type interactions 
with a "background gas" of other spin particles. These 
processes can lead to both local spin exchange and 
de-phasing. Here we model these processes again by 
a Lindblad- type master equation [3 , but we mention 
that a numerical treatment including memory effects in 



the environment, which gives rise to non-Markovian de- 
coherence, can be given [2 . 

We calculate the entanglement that is generated during 
the molecular motion, using the two-qubit measure of 
concurrence 

Master equation for the oscillating molecule in 
contact with bosonic heat baths.— The derivation of 
the master equation follows standard arguments used in 
reservoir theory (see e.g. [l]), but with some important 
modifications due to the time dependence of the problem. 

One should point out from the outset that the deriva- 
tion of the master equation rests on a series of assump- 
tions (including e.g. the Born and Markov approximation 
and the secular or rotating- wave approximation), neither 
of which we expect to be very well satisfied in real bio- 
logical systems. 

What the master equation does provide^ however, is 
a dynamical process that exhibits the essential features 
that we expect to be most relevant in our system of con- 
sideration: A process that is disentangling and that leads 
to de-coherence and thermalization in the subsystem, the 
two-spin molecule. Whether or not these processes fol- 
low, e.g., an exponential decay is not so essential for the 
main argument. 

The total Hamiltonian of the system and the environ- 
ment can be written in the form 

Htot {t) = Hm {t) ^Hb^ Hmb (5) 

where HM{t) is the time-dependent Hamiltonian (LI) of 
the two-spin molecule and Hb, Hmb describe the oscilla- 
tor bath and the molecule-bath interactions, respectively. 
We assume dissipative coupling, in which case the latter 
can be written in the form 

ffMB = ^<T("^®S„ (6) 

a 

where = B^^ denote the collective bath degrees of 
freedom that couple to the ath spin. The interaction 
constants have been absorbed in the ^^^s. 

Within the usual Born-Markov approximation, one ob- 
tains an equation of motion for the two-spin molecule 
which, in the interaction picture, has the form 

^p(t) = - J dsTrB[HMBit),[HMB(t-s),p{t)^pB]l 

(7) 

where pit) — TrBpMB{t) represents the reduced density 
matrix of the two-qubit molecule after tracing out the 
degrees of freedom the thermal bathes. 

To perform the secular or rotating wave approxima- 
tion, we expand the spin operators ai^^ in (p| into the 
basis of instantaneous eigenstates \ei{t)) with eigenvalues 
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ei{t) {i = 0, .., 3) of the system Hamiltonian HM{t), i.e. 

= E E slfit)\em{eM- (9) 

In the first line, ttie operators S^^^^^ describe intra- 
molecular transitions with frequency uj{t) and the sum- 
mation runs over all resonant transition frequencies; in 
the second line, S^^f (t) = (ei(t)|ai"^ |e,(t)) are the cor- 
responding transition matrix elements and the summa- 
tion runs over all states with matching energy eigenvalues 
ei(t), ej{t). In the interaction picture, which is used in 
the derivation of the master equation, we then write 



r(«) 



it) 



■■^uHt)sl%u{t) 



(10) 



where U (t) is the unitary evolution generated by the sys- 
tem Hamiltonian, U{t) = T ^-"^ So HM{s)ds ^ j (^e^otes 
the time ordering operator. 

In a situation where the system's evolution is adiabatic^ 
i.e. slow enough to avoid energy-changing transitions, the 
eigenstates will only pick up a dynamic phase (under the 
coherent evolution of the time-dependent system Hamil- 
tonian) 0. The time dependence of the spin operators 
in the interaction picture then acquires the simple form 



it) 



Upon inserting this into ^ and applying (a general- 
ization of) the secular approximation [1 , we obtain a 
master equation of the form (L2) with implicitly time- 
dependent Lindblad operators [5 . The properties of the 
baths thereby enter through the following Fourier-type 
transform of the bath correlation functions 

poo ^ 

r^0H-),t]= dse^/-.-W^«(i?t(,)B^(0)) (12) 



which is different from the exact Fourier transform ob- 
tained in case of a time-independent system Hamiltonian. 
For independent baths, we only get a contribution for 
a = (3. The Markov approximation assumes that the 
correlation functions {Bl^{s)Ba{0)) decay fast compared 
to the relaxation time, which means that their real and 
imaginary part can essentially be replaced by the delta 
function d{t) and its time derivative S\t), respectively. 
It is therefore consistent to apply, for small values of 5, 



the following approximation for the integral in (12) 



j: 



uj{t')dt' ^ uj{t)s 



In summary, we can thus write 



it) 



(13) 



(14) 



where the function r^(^) now depends only on the value 
of the frequency u{t) at the time t. Upon transform- 
ing back to the Schrodinger picture, \ei{0)) in the tran- 
sition operators cri^\t) are mapped to |ei(t)), and we 
finally obtain a master equation with the implicitly time- 
dependent Lindblad generators 



slf{t)\eM){e,{0)\- 



(11) 



= La(w(t)) 



a/2 



Ay(t)=w(t) 



(15) 

The index /i = {a,a;(t)} in Eq. (L2) of the Letter runs 
here over a = 1,2 and over the allowed values of u;{t). 
The relevant quantity of the heat bath which enters F^^^) 
is its spectral density function. If we assume an Ohmic 
spectral density with infinite cut-off frequency, we obtain 



F,(,) =/^Cj(t)(l+7V,(,),;3) 



(16) 



where N^(t),f3 is the bosonic distribution function at in- 
verse temperature /3, i.e., N^(t),f3 = l/{e^^^'^^ — 1). The 
master equation (L2) with these Lindblad operators was 
used to calculate the data shown in Fig. L3. 

Concurrence and static thermal entanglement.— 

To measure the dynamic entanglement generated dur- 
ing molecular oscillation, we compare it with the ther- 
mal equilibrium state when the molecular configuration 
is fixed at any distance. In such a "static configuration" , 
both the Hamiltonian and the Lindblad operators are 
time independent. Furthermore, the derived quantum 
master equation is then mixing and the molecule will al- 
ways be driven to its thermal equilibrium state [1 corre- 
sponding to the reservoir temperature /3 = 1/T. For each 
specific molecular configuration, with fixed spin-spin in- 
teraction strength J and local electric or magnetic fields 
the thermal equilibrium state reads 



Pth 



1 

1 



/ roo ro3 \ 

rii ri2 

r2i r22 

V ^30 rss J 



(17) 



where Z = Tr{e~^^^) is the partition function and the 
matrix representation refers to the computational prod- 
uct basis. The non-zero entries of the above matrix are 
given by 



^00 



2sinh {SP)/{l^r]^), 



(18) 

rii = r22 = cosh (J/3), (19) 
rss = e-^^ + 2sinh(^/3)/(l + 7?2), (20) 
ro3 = -Jsmh{S(3)/S, ri2 = - sinh ( J/3) (21) 



where Z = 2[cosh {S(3) + cosh (J/3)], £ = {4B^ + J^) 
and T] = {S — 2B)/J. To quantify the two-qubit entan- 
glement, there exist various kinds of entanglement mea- 
surements. We choose the concurrence C [4 , which is 
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defined as C = max{0, Ai — A2 — A3 — A4}, where the 
A^s are the square roots of the eigenvalues of pp in de- 
creasing order [4 , with p = {ay (g) (jy)p*{(jy ^ ay). For 
the thermal equilibrium state (17), one obtains 



C{pth) 



- max{0, |ri2| - (roorsa)^/^, \ro3 \ - {rur22)^^^}. 

(22) 

Using the explicit expressions for ri2,roo,^33 one finds 
that |ri2| - (rooras)^/^ < and |ro3| - (rnrs)^/^ = 
^sinh {£(3) - cosh {J p). 
The static thermal entanglement can thus be written 

as 

2 I 71 

C{pth) = -max{0,y sinh(f/3)-cosh(J/?)}. (23) 

In order to illustrate how the static entanglement 
changes as the temperature increases, we calculate the 
first derivative of C with respect to p. After some 
straightforward calculations, it can be seen that 



dC{pth)/df3 > 



(24) 



which means that the static entanglement always de- 
creases as the temperature increases. This is consistent 
with our intuition that there exists a critical temperature 
Tc, above which no static entanglement can survive. In 
other words, for any fixed molecular configuration, en- 
tanglement will eventually vanish when the reservoir is 
too hot. In our Letter, we consider exactly such a sit- 
uation. The temperature of the environment is so high 
that the thermal state is separable (non-entangled) at all 
possible molecular configurations, i.e. T > max{Tc}. 

Spin-gas model for the environment.— The envi- 
ronment of bio-molecular systems is rather complex [6] 
and not yet fully understood, and some of its features 
may not be well-described by a thermal bath model of 
harmonic oscillators. To check whether the observed ef- 
fects are robust, we have also considered an alternative 
model - the so-called spin gas model - for the en- 

vironment. In this model, we assume collisions between 
the molecular spins and other, randomly moving spin- 
particles that constitute the environment. The collisions 
induce local energy dissipation (i.e. spin exchange) and 
de-phasing, which leads to de-coherence and, if left alone, 
quickly destroys all entanglement between the molecular 
spins. [12 

The spin-gas model of the environment has been de- 
scribed in Refs. 00. For Ising-type interactions, one can 
calculate the time evolution of small subsystems exactly, 
for environments consisting of up to 10^ — 10^ particles, 
without any approximations. Under certain conditions, 
one can again derive a master equation by considering 
the effect of random collisions on a coarse-grained time 
scale [3]. For the present purpose, we will employ such 
a phenomenological description, but we emphasize that 




Figure 5: Ground state population po (left) and dynamic 
entanglement C (right) vs. the molecular configuration d. 
Upper: Bosonic thermal bath at temperature T = 1 and 
K = 0.01; Lower: Spin gas model with s = 0.16 and 7 = 0.025. 
The oscillation parameters are the same as in Fig. L3 of our 
Letter. The first cycle and the asymptotic cycle are indicated 
by the blue and red curves respectively. The green curves 
(left) represent the instantaneous steady state corresponding 
to each molecular configuration. The blue dot and arrows 
indicate the starting point and the evolution direction. 



non-markovian and collective effects in the environment 
can be taken into account [2 . 

The effect of random collisions leads to Lindblad gen- 
erators of the form Li G {Lg^\L^^^\a = 1,2} with 



where = (^r ± if^r )/2 

are the Pauli ladder oper- 
ators for a two-level system. The resulting master equa- 
tion |3] describes local energy gain and loss processes 
("spin exchange") with the effective rate 7 > 0, while s 
is related to the temperature and determines the equi- 
librium distribution of the local excitation 0. Without 
loss of generality, we may assume that the 7 > and 
0<5<l/2. If 5 is larger than a critical value 5c, no 
static entanglement can exist [5^, similar as in the case 
of the bosonic heat bath in the previous section. 

Even though this model is quite different from the 
bosonic heat bath, it has certain features in common, for 
example, it is disentangling and mixing. Remarkably, we 
find the same phenomenon as in case of the bosonic heat 
bath, namely a persistent recurrence of fresh entangle- 
ment in a regime where the static entanglement vanishes 
for all molecular configurations. 

In Fig. [5] we compare the evolution of the oscillating 
molecule for the spin gas model and the bosonic heat 
bath model. The upper panel of Fig. [5] reproduces the 
left 3D plot in Fig. L3 (left) of our Letter by projecting 
it onto two dimensional curves. The lower panel shows 



(a) 



and L 



(a) 
d 



^7(1 - s)a. 



(a) 



(25) 
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the same evolution for the spin gas. It can be seen that 
the qualitative features are robust: During the first oscil- 
lation, entanglement builds up when the spins approach 
each other, while the evolution subsequently converges 
towards an asymptotic cycle on which the entanglement 
periodically recurs. This observation strengthens our 
claim that this feature is robust and does not seem to 
depend on the detailed features of the environment. 

Competing effects of environmental noise on 
dynamic entanglement.— Another benefit of the phe- 
nomenological spin-gas model is that it allows one to en- 
ter the regime of short oscillation periods and to clearly 
illustrate the competition between the constructive and 
the destructive effects of the environmental noise. [13] In 
Fig. L4, we have plotted the maximal value of entangle- 
ment that is assumed on the asymptotic cycle, in the case 
where s = 0.2 > Sc (i.e. the equilibrium state is sepa- 
rable for all possible molecular configurations). The left 
plot displays the maximally achievable entanglement for 
different oscillation periods. It can be seen that the oc- 
currence of dynamic entanglement strongly depends on 
the oscillation period; there are competing effects of the 
environmental noise which give rise to an optimal oscil- 
lation period where the effect dynamic entanglement is 
most pronounced. 




10 15 20 30 50 ■ 10 15 20 30 50 



d d 

Figure 6: Ground state population po (left) and dynamic en- 
tanglement C (right) vs. the molecular configuration d (loga- 
rithmic coordinate). The oscillation periods are r = 6 (green), 
20 (red) and 100 (blue) respectively. The dot and arrow in- 
dicate the starting point and the evolution direction. Right: 
The green curve (with positive values) only represents the 
first cycle, while subsequent cycles collapse onto the abscissa, 
without entanglement. The other parameters are the same as 
in Fig. L4 (left) of our Letter. 

To understand the results in Figs. L4 and [6] better, 
it is worth looking at the detailed time evolution of the 
entanglement and the ground state population for three 
typical oscillation periods. 

First, consider a very long oscillation period, e.g. 
r = 100. Here, the molecule is almost completely re- 
set by thermalization to equilibrium when the spins are 
spatially separated (distant configuration), with a large 
ground state population of up to ~ 65%, see Fig.[6](bl ue 
curve). Since in this regime the coherent evolution of 
the molecule is adiabatic, the population of the instan- 



taneous eigenstates of the system Hamiltonian remains 
approximately constant, while the off-diagonal elements 
remain negligible. When two spins come closer, they 
start interacting and the ground state becomes entan- 
gled (entanglement generation regime). If there were no 
dissipation, the high population of the ground state alone 
would be sufficient to generate entanglement. However, 
while the spins approach each other, the energy sepa- 
ration between the lowest lying levels decreases and the 
dissipation starts re-populating the levels. This drives 
the molecular state into the separable regime and dimin- 
ishes its entanglement, with only little entanglement sur- 
viving. For a moderate oscillation period, e.g. r = 20 
(red curve), the dissipation still has enough time to reset 
the system while it passes through the distant configu- 
ration, with a ground state population similar as for as 
r = 100. In the entanglement generation regime, how- 
ever, the destructive effect of the dissipation is now much 
smaller than for long oscillation period, which leads over- 
all to more entanglement. Finally, for a very short oscil- 
lation period, e.g. r = 6, the destructive effect during 
the entanglement generation regime is even smaller. On 
the other hand, the reset effect in the distant configu- 
ration is greatly suppressed, since the system does not 
have enough time to thermalize and to repopulate the 
ground state. Thus, even though the transient entan- 
glement is larger in the first period, as expected, it will 
diminish in subsequent runs and cannot be sustained on 
the asymptotic cycle, due to the lack of an effective reset 
mechanism. 

Generic features of persistent dynamic entan- 
glement in noisy environment. — In the main text, 
we have used a simple model, namely a harmonically os- 
cillating molecule in an Ohmic thermal bath, to demon- 
strate the essential mechanism for persistent dynami- 
cal entanglement to occur in a de-coherent environment 
where no static entanglement can exist. Here we show 
that neither the harmonic oscillatory motion nor the 
Ohmic thermal bath are indispensable for such an effect. 
The existence of persistent dynamic entanglement is to 
a very large extent independent of the precise details of 
the classical motion and thermal environment. 

In Fig. [g we consider a model where the spins move 
towards (away) from each other with a constant speed, 
and observe similar results as for the harmonic oscillatory 
motion. The same effect can also be seen in a scenario 
of stochastic movements [5]. In short, the detailed char- 
acteristics of the entanglement (how much entanglement, 
how does it vary with time, etc.) depend on the driving 
oscillation, but the very existence of persistent entangle- 
ment is generic. All that is needed is that the classical 
motion obeys two conditions: (i) is adiabatically slow 
but sufficiently fast compared to de-coherence and (ii) it 
spends long enough time at the far end for thermalization 
to occur. 
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Figure 7: Ground state population pg and entanglement C 
vs. the molecular configuration characterized by spin-spin 
distance d for the bosonic heat bath. The spins move towards 
(away) from each other with a constant speed. The other 
parameters are the same as Fig. L3. 




Figure 8: Dynamic entanglement C vs. the molecular config- 
uration d for the sub-ohmic (left, s = 0.8) and supra-ohmic 
(right, s = 1.2) thermal bath. The other parameters are the 
same as Fig. L3. The blue dot and red arrows indicate the 
starting point and the evolution direction. The first cycle and 
the asymptotic cycle are indicated by the blue and red curves, 
respectively. 

Regards the model for the bath, we have so-far used the 
Ohmic bath (as well as the spin-gas model) as a exam- 
ple. However, our results are also valid for other forms 
of spectral densities. It can be seen from Fig. [8] that 
the effect of persistent dynamic entanglement is not re- 
stricted to the simple Ohmic bath, but occurs also for 
the sub-Ohmic and supra-Ohmic bath, i.e. for a spectral 
density ~ uj^ with 5 < 1 or s > 1. We have also found 
that the present mechanism works even with the spec- 



tral density from the solvent and protein environment 
|8], e.g. a spectral density in the form of ~ ^. Finally, 
by using the numerical method of quasi-adiabatic propa- 
gator path integral [9l |10] , we have extended the results 
to the non-Markovian environment with finite memory 
time, and still see the generic effect discussed in the main 
text. More details will be presented in [5]. 
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